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The formulation of gauge theories on compact Riemannian manifolds with 
boundary leads to partial differential operators with Gilkey-Smith boundary con- 
ditions, whose peculiar property is the occurrence of both normal and tangential 
derivatives on the boundary. Unlike the standard Dirichlet or Neumann boundary 
conditions, this boundary-value problem is not automatically elliptic but becomes 
elliptic under certain conditions on the boundary operator. We study the Gilkey- 
Smith boundary-value problem for Laplace-type operators and find a simple crite- 
rion of ellipticity. The first non-trivial coefficient of the asymptotic expansion of 
the trace of the heat kernel is computed and the local leading asymptotics of the 
heat-kernel diagonal is also obtained. It is shown that, in the non-elliptic case, the 
heat-kernel diagonal is non-integrable near the boundary, which reflects the fact 
that the heat kernel is not of trace class. We apply this analysis to general linear 
bosonic gauge theories and find an explicit condition of ellipticity. 



1 0n leave of absence from Research Institute for Physics, Rostov State University, Stachki 194, 344104 
Rostov-on-Don, Russia. 



I. G. Avramidi and G. Esposito: Heat kernel asymptotics 



1 



1 Introduction 

Elliptic differential operators on manifolds have proved to play a crucial role in mathe- 
matical physics. In particular, the main objects of interest in quantum field theory and 
statistical physics, such as the effective action and the partition function, are described 
by the functional determinants, or, which is equivalent, by the zeta-function and the 
heat kernel of self-adjoint elliptic differential operators. Of particular importance are, 
of course, the operators of Laplace type or Dirac type. |] In the case of manifolds with 
boundary, one has to impose some boundary conditions in order to make a (formally self- 
adjoint) differential operator self-adjoint and elliptic. Indeed, the boundary conditions 
are additional ingredients in the theory which have not been fixed a priori, and the choice 
of boundary conditions is, by no means, unique. There are many admissible boundary 
conditions that guarantee the self-adjointness and ellipticity of the problem. The simplest 
boundary conditions are the classical Dirichlet and the Neumann ones. In the Dirichlet 
case one sets to zero at the boundary the value of the field, whereas in the Neumann case 
the normal derivative of the field is set to zero at the boundary. There exist also slight 
modifications of the Neumann boundary conditions (called Robin boundary conditions in 
physical literature) when the normal derivative of the field at the boundary is not set to 
zero but is proportional to the value of the field at the boundary. An even more general 
scheme, called mixed boundary conditions, applies to the operators acting on sections 
of some vector bundles. It is then possible to mix the Dirichlet and Robin boundary 
conditions by using some projectors, i.e. a part of the field components satisfy Dirichlet 
boundary conditions and the remaining part satisfies Robin boundary conditions. 

However, this is not the most general scheme. A much more general setup for the 
boundary- value problem was developed by Gilkey and Smith. [EJ] They put forward some 
boundary conditions that are still local but include both normal and tangential derivatives 
of the fields at the boundary. In this paper we are going to study the Gilkey-Smith 
boundary-value problem for operators of Laplace type. Unlike the Dirichlet or Neumann 
boundary-value problems, such a boundary-value problem is not automatically elliptic. 
Therefore, following Refs. 1,3, we find, first, a criterion of (strong) ellipticity, which 
provides an explicit simple condition on the boundary operator. Then we construct the 
parametrix to the heat equation in the leading approximation and compute the first 
non-trivial (next to leading) term Ai/ 2 in the asymptotic expansion of the trace of the 
heat kernel. We also discuss what happens when the boundary-value problem is not 
strongly elliptic. Last, we study the problem of ellipticity in linearized gauge theories on 
manifolds with boundary. The attempt to preserve gauge invariance on manifolds with 
boundary fixes the boundary conditions and leads exactly to a Gilkey-Smith boundary- 
value problem. As is shown in Ref. 1, while Yang-Mills as well as Rarita-Schwinger 
theories are automatically elliptic, quantum gravity based on the Einstein action turns 
out to be not elliptic, if the Gilkey-Smith boundary-value problem is studied. 
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2 Gilkey— Smith Boundary- Value Problem 

Let (M, g) be a smooth compact Riemannian manifold of dimension m with smooth 
boundary dM. Let g be the positive-definite Riemannian metric on M and g be the 
induced metric on dM. By using the inward geodesic flow, we identify a narrow neigh- 
bourhood of the boundary dM with a part of dM x R + and define a split of the cotangent 
bundle T*(M) = T*(dM) © T*(R). Let z = (£*), with z = 1, 2, . . . , m - 1, be the lo- 
cal coordinates on dM and r be the normal geodesic distance to the boundary, so that 
N = d r = d/dr is the inward-pointing unit normal vector field to the boundary. Near 
dM we choose the local coordinates x = = (ct, r), with \x = 1,2, ... ,m, and the split 
f = (Q = ((,uj) G T*{M), where ( = (£,) e T*(dM) and w 6 R. With our notation, 
Greek indices range from 1 through m and lower case Latin indices range from 1 through 
m — 1. 

Let V be a (smooth) vector bundle over the manifold M and C°°(V, M) be the space 
of smooth sections of the bundle V. Let V* be the dual vector bundle and E : V — > V* 
be a Hermitian non-degenerate metric, E* = i£, that determines the Hermitian fibre 
scalar product in V. Using the invariant Riemannian volume element dvol(x) on M one 
defines a natural L 2 inner product (, ) in C°°(V,M), and the Hilbert space L 2 (V,M) as 
the completion of C°°(V, M) in this norm. 

Let V T *^ M - ) be the Levi-Civita connection on M and V v be the connection on the 
vector bundle V compatible with the metric E. Then we define, as usual, \7 T *( M )® V — 
yT*(Af) (g) ]__)_]_ y^. Moreover, we will often denote just by V the total covariant 
derivative without mentioning the bundle it is acting on. The notation V will be used to 
denote the covariant tangential derivative on the boundary. 

Let, further, tr g = g®l be the contraction of sections of the bundle T* \M)®T* (M)®V 
with the metric on the cotangent bundle T*(M), and Q G C°°(End (V), M) be a smooth 
self-adjoint endomorphism of the bundle V, i.e. Q = E~ X Q}E = Q. Then a Laplace-type 
operator F : C°°(V, M) — > C°°(V, M) is a second-order differential operator defined by 

F = _tr g V T * (A/) ® y V y + Q. (2.1) 

Let us define the boundary data by 

where ^(v 9 ) = V^I^m and - ?/'i((^) = ^NfldM are the restrictions to the boundary of the 
sections if G C°°(V, M) and their normal derivatives. Let the vector bundle W over 
dM be the bundle of the boundary data. W consists of two copies of the restriction 
of V to dM and inherits a natural grading^ W = Wo © W\, where Wj represents 
normal derivatives of order j, and, therefore, dimp^ = 2 dim V. The bundles Wo and W\ 
have the same structure, and hence in the following they will be often identified. Let 
W' = Wq © W[ be an auxiliary graded vector bundle over dM such that dim W' = dim V. 
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Let B : C°°(W,dM) -> C°°(W',dM) be a tangential differential operator on dM. The 
boundary conditions then read 

Bip((p) = 0. (2.3) 

Now let II be a self-adjoint projector acting on W . Since the bundle W\ is identified 
with W the projector n acts also on W x . Let Y e C°°(T(dM) ® End (W ),dM) be an 
anti-self-adjoint endomorphism- valued vector field on dM orthogonal to II, i.e. P = — T l , 
HP = r*II = 0, and S £ C°°(End (Wo)j 9M) be a self-adjoint endomorphism orthogonal 
to II, i.e. S = S, IIS = SU = 0. Using these objects we define a first-order self-adjoint 
tangential differential operator A : C°°(W ,dM) -> (7°°(Wo,9M) by 

a = (/ - n) {-(rv* + v t r) + 5} (/ - n), (2.4) 

that is obviously orthogonal to II: LTA = ALT = 0. Hereafter / is the identity endomor- 
phism of the vector bundle V. The Gilkey-Smith boundary operator is expressed in terms 
of these geometric objects by 

n o 



Ma i-a)- < 2 - 5 > 

It is not difficult to see that the Gilkey-Smith boundary-value problem incorporates all 
standard types of boundary conditions. Indeed, by choosing II = / and A = T = S = 
one gets the standard Dirichlet boundary conditions, by choosing 11 = 0, T = 0, A = S = I 
one gets the standard Neumann boundary conditions. More generally, the choice T = 
and A = S corresponds to the mixed boundary conditions mentioned in Sec. 1. 



3 Strong Ellipticity 

Integration by parts shows [|I|] that the Laplace- type operator F given in Eq. (EOT) en- 
dowed with the Gilkey-Smith boundary conditions ( |2.3j ) is symmetric, meaning that 
(<Pi,F(p 2 ) = (F(pi,ip 2 ) for any two smooth sections y?i,y?2 £ C°°(V,M) satisfying the 
boundary conditions Bip((pi) = Btp({p 2 ) = 0. However, it is not automatically elliptic. 
Now we are going to determine under which conditions the Gilkey-Smith boundary- value 
problem for a Laplace-type operator is strongly elliptic. 0] 

First of all, the leading symbol of the operator F should be elliptic in the interior of 
M. Let, hereafter, A be a complex number which does not lie on the positive real axis, 
A G C — R + (R + being the set of positive numbers). Then by using the leading symbol 
of the operator F, i.e. ctl(F;x, £) = |£| 2 • /, with |£| 2 = g^ix)^^, we find easily 

det (a L (F; x, - A • /) = (|£| 2 - A) dimy ? for (£, A) + (0, 0). (3.6) 

Thus, the leading symbol of the operator F is elliptic. 

Second, the so-called strong ellipticity condition should be satisfied. |2], £§] As we already 
noted above, there is a natural grading in the vector bundles W and W which reflects 
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simply the number of normal derivatives of a section of the bundle. The boundary 
operator B in Eq. ( |2.5|) is said to have the graded order 0. Its graded leading symbol is 
defined by[@, 1 

"•<*'> s (a? ( (/-n))> < 3 - 7 > 

where r • ( = T^Q. To define the strong ellipticity condition we take the leading symbol 
ctl(F] x, r, (, lj) of the operator F, substitute r = and uj — > — id r and consider the 
following ordinary differential equation for a section (p G C°°(V, dM x R+): 

[a L (F; x, 0, C, -«9 r ) - A • J] p(r) = 0, (3.8) 

with an asymptotic condition 

lim <p(r) = 0. (3.9) 

1 >CX) 

The boundary- value problem (F, B) is said to be strongly elliptic^ Q with respect to 
the cone C — R + if for every pair ((, A) ^ (0,0), and any ip' G C 0O (W, 9M) there is a 
unique solution <p of the equation (|3.8|) satisfying the asymptotic condition (|3.9|) and the 
boundary condition 

(r g (B F ){x,CM<p) = i!/, (3.10) 

where 4>((p) G C°°(W, dM) are the boundary data defined by (|2.2|) . 

For a Laplace-type operator this definition leads to the following theorem. [JTJ] 

Theorem 1 The Gilkey-Smith boundary-value problem (F, B) is strongly elliptic with 
respect to C — R + if and only if the matrix \(\I — iT-( is positive-definite, i.e. — iT ■( > 
0, for any non-vanishing (. A sufficient condition for strong ellipticity is: 

|C| 2 / + (r-o 2 >o. (3.ii) 



4 Asymptotic Expansion of the Heat Kernel 

For t > the heat semi-group exp(— tF) : L 2 (V,M) —>■ L 2 (V,M) of the strongly elliptic 
boundary-value problem {F,B) is well defined. The kernel of this operator, U(t\x,y), 
called heat kernel, is a section of the tensor product of the vector bundles V and V* over 
the tensor-product manifold M x M, defined by the equation 

(dt + F)U(t\x,y) = (4.12) 

with initial condition 

U(0 + \x,y) = 5(x,y), (4.13) 

where 5(x, y) is the covariant Dirac distribution. Moreover, a boundary condition is 
imposed, i.e. 

BiP[U(t\x,y)] = 0, (4.14) 
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and the symmetry condition holds 

U(t\x,y) = U(t\y,x). (4.15) 

Hereafter all differential operators as well as the boundary data map act on the first 
argument of the heat kernel, unless otherwise stated. 

It is well known|4]] that the heat kernel U(t\x,y) is a smooth function near diagonal 
of M x M and has a well defined diagonal value U(t\x, x), and that the L 2 trace 

Tr L 2 exp(-tF) = f dvol{x)tr v U{t\x,x), (4.16) 

JM 

has an asymptotic expansion as t — > + 

Tr L2 exp(-tF) ~ (Ant)^ 2 £ t k ' 2 A kl2 (F, B). (4.17) 

fc>0 

Here Ak/ 2 {F,B) are the famous global heat-kernel coefficients (sometimes called also 
Minakshisundaram-Plejel or Seeley coefficients). The zeroth-order coefficient is very well 
known: 

A Q = J dvol(x)tr y J = vol (M) • dim(V). (4.18) 

M 

It is independent of the operator F and of the boundary conditions B. The higher order 
coefficients have the following general form: 

A k/2 (F,B) = J dvol{x)tr v a k/2 {F\x) + J dvol(x)tr V b k/2 {F, B\x), (4.19) 

M dM 

where a k / 2 {F\x) and b k /2{F, B\x) are the (local) interior and boundary heat-kernel coef- 
ficients. The interior coefficients a k /2{F\x) do not depend on the boundary conditions. 
Moreover, it is well known that they vanish for half- integer order, a k+ i/ 2 = 0.0] The in- 
teger order coefficients a k (F\x) are calculated for Laplace-type operators up to 04. The 
boundary coefficients bk/ 2 (F, B\x) do depend on both the operator F and the boundary 
operator B. They are far more complicated because in addition to the geometry of the 
manifold M they depend essentially on the geometry of the boundary dM. For Laplace- 
type operators they are known for the usual boundary conditions (Dirichlet, Neumann, 
or mixed version of them) up to 65/2 -M, For Gilkey-Smith boundary conditions only 
some special cases have been studied in the literature. || ||, |TT], [12| In this paper we 
evaluate the next-to-leading coefficient Ai/ 2 (F, B), following our recent work.[|T] 



5 Parametrix: General Setup 

In this section we show how one can construct an approximation to the heat kernel 
U(t\x, y) for t — > + near the diagonal, i.e. for x close to y. First of all, we decompose 
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the heat kernel into two parts 

U(t\x, y) = U^x, y) + U B (t\x, y). (5.20) 

Then we construct different approximations for Uoo and Ub- The first part U^^x, y) is 
approximated by the usual asymptotic expansion of the heat kernel in the case of compact 
manifolds without boundary when x — > y and t — > + . This means that effectively one 
introduces a small expansion parameter e reflecting the fact that the points x and y are 
close to each other and the parameter t is small. This can be done by fixing a point x', 
choosing the normal coordinates at this point (with g^pix') = 5^ u ) and scaling 

x — > x' + e(x — x'), y — > x' + e(y — x'), t — > e 2 t, (5-21) 

and expanding into an asymptotic series in e. This construction is, by now, quite 
standard^] and we do not repeat it here. One can also use a completely covariant 
method, || which leads to the result 

~ (4vrt)- m / 2 exp (-^&) J2t k a k (x,y), (5.22) 



fc>0 



where d(x,y) is the geodesic distance between x and y and a k (x,y) are the off-diagonal 
heat-kernel coefficients. These coefficients satisfy certain differential recursion relations 
which can be solved in form of a covariant Taylor series near diagonal. || On the diagonal 
the asymptotic expansion of the heat kernel reads 

£/oo(t|z,aO ~ (A7it)- m ^J2t k a k (F\x), (5.23) 

fc>0 

where a k (F\x) = ak(x,x). As we noted above, the explicit formulae for the diagonal 
values of a k are known up to k = 4. 0] This asymptotic expansion can be integrated over 
the manifold M to get 

J dvol(a;)tr yU^x, x) ~ (Ant)-^ 2 ^ t k J dvol(x)tr v a k {F\x). (5.24) 

M M 

Thus, integrating the diagonal of gives the interior terms in the heat-kernel asymp- 
totics fl£ia|) . 



For a strongly elliptic boundary- value problem the diagonal of the boundary part 
UB{t\x,x) is exponentially small as t — > + if x dM, i.e. of order ~ exp(— r 2 (x)/t), 
where r(x) is the normal geodesic distance from x to the boundary. Thus, it does not 
contribute to the asymptotic expansion of the heat-kernel diagonal outside the boundary 
as t — > + . This implies that the asymptotic expansion of the total heat-kernel diagonal 
outside the boundary is determined only by 

U(t\x, x) ~ (47rf)~ m/2 t k a k {F\x), x £ dM. (5.25) 

k>0 
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The coefficients of the asymptotic expansion as t — > + of the diagonal of the boundary 
part UB{t\x,x) behave near the boundary like the one-dimensional Dirac distribution 
S(r(x)) and its derivatives. Thus, the integral over the manifold M of the boundary 
part UB{t\x,x) has an asymptotic expansion as t — > + with non-vanishing coefficients 
in form of integrals over the boundary. The local boundary coefficients bk/2 contribute, 
after integration over the boundary, to the global heat-kernel coefficients, according to 
Eq. ( }4.19|) . It is well known that the coefficient A1/2 is a purely boundary contribution. 0] 



It is almost obvious that it can be evaluated by integrating the fibre trace of the boundary 
contribution Ub of the heat kernel to leading order. 

Of course, Uoo is obtained without taking into account the boundary conditions. There- 
fore, it satisfies approximately the equation ( 4.12|) but does not satisfy the boundary con- 



ditions (|4.14p . This implies that the compensating term UB(t\x,y) should be defined by 
the equation 

(d t + F)U B (t\x,y) = (5.26) 



with the initial condition 
and the boundary condition 



U B (0 + \x,y)=0, (5.27) 



Bi, [U^tlx, y) + U B (t\x, y)\ = 0. (5.28) 

The compensating term U B {t\x,y) is important only near the boundary where it behaves 
like a distribution when t — > + . Since the points x and y are close to the boundary the 
coordinates r(x) and r(y) are small separately, hence not only the difference [r(x) —r(y)] 
but also the sum [r(x) + r(y)] is small. This means that we must additionally scale 
r(x) — > er(x) and r(y) — > er(y). By contrast, the point x! is kept fixed on the boundary, 
so the coordinates x' do not scale at all: x' — > x' . 

Thus, we shall scale the coordinates x = (x,r(x)), y = (y,r(y)), and the parameter t 
according to 

x — > x 1 + e(x — x f ), y — > x' + e(y — x'), (5.29) 
r(x) — > er(x), r(y) — > er(y), t — > e 2 t. (5.30) 
The corresponding differential operators are scaled by 

1 „ 1 1 

d -> -d, d r -> ~d r , dt -> -rd t . (5.31) 

see 1 

We call this transformation just scaling and denote the scaled objects by an index 
e, e.g. U B . The scaling parameter e is considered as a small parameter in the theory 
and we use it to expand everything in power series in e. We do not take care about the 
convergence properties of these expansions and take them as formal power series. In fact, 
they are asymptotic expansions as e — > 0. At the very end of calculations we can set 
e — \. The non-scaled objects, i.e. those with e = 1, will not have the index s, e.g. 
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U B \ £= \ = Ub- Another way of doing this is by saying that we expand all quantities in the 
Taylor series in the boundary coordinates x and y about the point x' with the coefficients 
being homogeneous functions of r(x), r(y) and t. 

First of all, we expand the scaled operator F £ in power series in e 

F £ ~-£e n - 2 F n , (5.32) 

n>0 

where F n are second-order differential operators with homogeneous symbols. The bound- 
ary operator requires a more careful handling. Since half of the boundary data (|2.2j) 
contain normal derivatives, formally ipo — y?|r=o and ipi = d r (p\ r= o would be of different 
order in e. To make them of the same order we have to assume an additional factor e in 
all ipi € C°°(Wi, DM). Thus, we define the graded scaling of the boundary data map by 

so that the boundary data map if) does not scale at all. This leads to an additional factor 
e in the operator A determining the boundary operator B of Eq. (|2.5| ). Thus, we define 
the graded scaling of the boundary operator by 

B '=(Z /-%•)• (5 - 34 > 

which has the following asymptotic expansion in e: 

B?~Y,e n B {n) , (5.35) 

n>0 

where Bt n ) are first-order tangential operators with homogeneous symbols. At zeroth 
order we have 

F = - d\ (5.36) 

where 

d 2 = g jk {x')djd k , A = T j (x')d j , U = U(x'). (5.38) 

Note that all leading-order operators Fq, Biq\ and Ao have constant coefficients and, 
therefore, are very easy to handle. This procedure is called sometimes "freezing the 
coefficients of the differential operator" . 

The subsequent strategy is rather simple. Expand the scaled heat kernel in e and 
substitute into the scaled version of the equation (|5.26| ) and of the boundary condition 
([5.28 ). Then, by equating the terms of the same order in e one gets an infinite set 



of recursive equations which determine all Ubu)- The f/oo(n) are obtained simply by 
expanding the scaled version of ( |5.22j ) in power series in e. 
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6 Parametrix: Leading Order 

In this section we determine the parametrix of the heat equation to leading order, i.e. 
Uo = {7oo(o) + Ub(q)- As we already outlined above, we fix a point x' G dM on the 
boundary and the normal coordinates at this point (with §ik{x') = 5^), take the tangent 
space T(dM) and replace the manifold M by M = T(dM) x R + . By using the explicit 
form of the zeroth-order operators F , B and A given by ( |5.36| )-( |5.38| ) we obtain the 
equation 



and the boundary conditions 



d t -d 2 r -d 2 )U (t\x,y) = 0, 



U Q U (t\x,y) 



r(x)=0 

(/-IIo) (d r + iTid 3 )u Q (t\x,y) 



r{x)=0 



0. 



(6.39) 

(6.40) 
(6.41) 



where n = n(rr'), Tq = T^x'). Hereafter the differential operators always act on the first 
argument of a kernel. Moreover, for simplicity of notation, we will denote n and T just 
by n and T- 7 and omit the dependence of all geometric objects on x'. To leading order this 
cannot cause any misunderstanding. Furthermore, the heat kernel should be bounded, 



lim U (t\x,y) 



lim U {t\x,y) = 0, 

r(y)-*oo 



(6.42) 



and symmetric, 

U (t\x,y) = U (t\y,x). (6.43) 

To solve the above boundary-value problem we use the Laplace tarnsform in t and, 
since it has constant coefficients, the Fourier transform in (x — y). Therefore, it reduces 
to an ordinary differential equation of second order in r on R + , which can be easily solved 
taking into account the boundary conditions at r = and r — > oo. Omitting simple but 
lengthy calculations we obtain 



U (t\x,y) 



ui+ioo 



R» 



, (27T) 



m— 1 



dX 



e -iA+*(«-fi) G ( A | C) r(x)) r{y))j (6.44) 



2m 



where w is a negative constant and G is the leading-order resolvent kernel in momentum 
representation. It reads 



G(X\C,u,v) 



exp 



\u 



ICP 



A 



+ 



/ - 2n + 2ir • c /VICI 2 - a - zT • c 



exp 



u 



ICP 



(6.45) 
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where Re y |C| 2 — A > 0. Now, by scaling the integration variables A — > X/t and ( — > C/Vt 
and shifting the contour of integration over A (w — > w/t, which can be done because the 
integrand is analytic in the left half-plane of A) and using the homogeneity property of 
the resolvent kernel we obtain immediately 



U (t\x,y) = {Ant) 



-m/2 



w+ioo 



R» 



7r (m-l)/2 



G A 



c, 



exp { < 



r(x) r{y) 



Vt 



Vt 
ICI 2 



Vt 



(6.46) 



Next, let us change the variable A according to A = |£| + uj . In the upper half-plane, 
Imw > 0, this change of variables is single- valued and well defined. Under this change 
the cut in the complex plane A along the positive real axis from \(\ 2 to oo, i.e. ImA = 
0, \(\ 2 < Re A < oo, is mapped onto the whole real axis Imu; = 0, — oo < Reu; < +oo. 
The interval ImA = 0, < Re A < \(\ 2 on the real axis of A is mapped onto an interval 
Reu; = 0, < Imu; < \(\, on the positive imaginary axis of uj. As a function of uj the 
resolvent G is a meromorphic function in the upper half plane, Ima> > 0, with simple 
poles on the interval Reu; = 0, < Imw < on the imaginary axis. The contour of 
integration in the complex plane of uj is a hyperbola going from (e i3,r / 4 )oo through the 
point uj = \J\C\ 2 — w to (e J7r / 4 )oo. It can be deformed to a contour C that comes from 
— oo + is, encircles the point uj = i\(\ in the clockwise direction and goes to +oo + ie, 
where e is an infinitesimal positive parameter. The contour C does not cross the interval 
Recj = 0, < Imu; < on the imaginary axis and is above all the singularities of the 
resolvent G. 

After such a transformation we obtain 



U (t\x,y) = (47rf)- m / 2 



7r(m-l)/2 



exp {-\cr+K 



(x 



y) 



R" 



Vt 



X 



c 



ICI 2 + ^ 2 



c 



r(x) r(y) 



(6.47) 



Substituting here G given in Eq. ( |6.45|) and computing Gaussian integrals over uj and ( 
we obtain the "free" part 

|2\ 



Uoo(o)(t\x,y) = (Airt) 



-m/2 



exp 



F - V\ 
At 



(6.48) 



and the boundary part 
U B (p){t\x,y) = 



(47rt) 



exp 



\x — y\ 2 + [r(x) + r(y)Y 
At 



(/ - 2H) 



-n(t\x,y)\, 



(6.49) 
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where 



Q(t\x,y) 



Rm-1 



71 



x / — exp S^-cu 2 + icu ^)+j{y)} | r . f + r . c) -i_ (6 5Q) 

c 



Herefrom we obtain easily the diagonal value of the heat kernel: 



U {0) (t\x, x) = (4vrt)- m / 2 |/ + exp f-^J (/ - 211) + * . 



(6.51) 



where 



•W = - 2 I ^iw/ i |^ l<l "" 2+2 "" r 'C(-/ + r-C)- 1 . (6-52) 



du 



Rm-1 C 

This function can be expressed further as 







$(z) = -2e~ z 1-2— tf(z), 
oz 



where 



*(*) = / ^iwi°°^ ex p{-KI 2 -(p + ^ + 2^r.c}. 

R m-1 



(6.53) 



(6.54) 



It is not difficult to show that, as z — ■> oo, the functions \I/ and $ are exponentially small: 

1 



v ; 2z 



/-^(/ + r 2 ) + o(z- 4 ) 



2z 2 



$(z) ~ — e * 2 



-r 2 + o(z~ 



where T 2 = g^TT- 7 . For z = 0, by using the change ( — > — £, we obtain 



*(o) = ^ / ^o^expf-icP-fr-C) 2 }. 



(6.55) 
(6.56) 

(6.57) 



Rm-1 



Note that this integral converges only when the strong ellipticity condition \(\ 2 I+(T-() 2 > 
is satisfied. 
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7 A 1/2 (F,B) Coefficient 

Now we take the diagonal U^(t\x,x) given by (|6.51| ) and integrate over the manifold 
M. Because the boundary part Ub(o) is exponentially small as r(x) —>■ oo we can in fact 
integrate it only over a narrow strip near the boundary, when < r(x) < 5. The difference 
is asymptotically small as t — > + . Doing the change of variables z = r/\/i we reduce the 
integration to < z < 5/ \ft. We see that as t — > + we can integrate over z from to oo. 
The error is asymptotically small as t — > + and does not contribute to the asymptotic 
expansion of the trace of the heat kernel. 
Thus, we obtain 

Tr i2 exp(-tF) = / dvol(x)tr v U Q (t\x, x) + 0(r m/2+1 ) 
Jm 

= (4irt)- m / 2 {A + ViA 1/2 (F,B) + O(t)}, (7.58) 
where A is given by ( |4.18| ) and 

A 1/2 (F,B)= J dvol(x)tr y 6i/ 2 , (7.59) 

8M 

with 



6i /2 = -^(/ + 2n) + 2^(0) (7.60) 
Now, using (|6.53|) and ( |6.57|) and the fact that ^(oo) = we get easily 



b 1/2 = -^(I + 2U) + V^ J ^TWexpf-lCp-^-C) 2 }. (7.61) 

Rm-l 

Note again that this integral converges only when the strong ellipticity condition is satis- 
fied, i.e. |C| 2 /+(r-C) 2 >0. 

Further calculations of general nature, without knowing the algebraic properties of the 
matrices r j , seem to be impossible. One can, however, evaluate the integral in form of an 
expansion in the matrices P. The integral over ( becomes Gaussian, which enables one 
to obtain 

h/2 = ^[i ~ 211 + 2 £ J^-^ ■ ..g i2n _ li2 T^ ■ • -r-)|. (7.62) 

Since our main result (|7.61[) is rather complicated, we now consider two particular 
cases of physical relevance. 

I. The first non-trivial case is when the matrices T l form an Abelian algebra, i.e. 

[r\r j ] = o. (7.63) 
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One can then easily compute the integral ( |7.61| ) explicitly and obtain 

61/2 = ^ {-/ - 211 + 2(1 + r 2 )- 1 / 2 } . (7.64) 

In the case T = we recover the familiar result for mixed boundary conditions. [|], [f| 
In the case IT = 0, this coincides with the result of Ref. 8, where the authors consid- 
ered the particular case of commuting P matrices (without noting this explicitly). 

II. A very important case is when the operator A is a natural operator on the boundary. 
Since it is of first order it can be only the generalized Dirac operator. In this case 
the matrices F satisfy a Dirac-type condition 

PF + FP = 2 g ij -7— — rr 2 , (7.65) 



m 



which leads to [|TJ 



-(m-l)/2' 



6 1/2 = ^j-/-2n + 2(/ + I -- TT r 2 ] }. (7.66) 

Note that this differs substantially from the result of Ref. 8, and shows again that 
the result of Ref. 8 applies actually only to the completely Abelian case, when all 
matrices F commute. Note also that, in the most interesting applications (e.g. in 
quantum gravity), the matrices P do not commute. ||10|| The result ( |7.61| ), however, 



is valid in the most general case. A particular realization of the above situation is 
the "pure" Dirac case when P = —^(m — 1)(J — II), where k is a constant. In this 
case we have W\ 



h/2 = ^-{-n+(/-n) [2(1 - K)- {m - 1]/2 - l] } . (7.67) 

Thus, a singularity is found at k — 1. This happens because, for k — 1, the strong 
ellipticity condition is violated (see also Ref. 9). Indeed, the strong ellipticity condition 

(r • o 2 + \C\ 2 i = ICI 2 [n + (i - k)(i - n)] > o, (7.68) 

implies in this case k < 1 (cf. Ref. 9). This is a general feature of the Gilkey-Smith 
boundary-value problem: the heat-kernel coefficients have singularities when the strong 
ellipticity condition is violated. H 
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8 Boundary Singularities in the Non-Elliptic Case 



Let us now consider the special case when V is a spin-tensor bundle and the boundary 
operator A is a natural operator, i.e. the matrices T- 7 form a representation of Spin(m). In 
other words, the matrices T- 7 can be constructed only from natural objects like the metric, 
the normal, the Dirac matrices and the frame. Then the eigenvalues of the leading symbol 
of the operator A, i.e. the matrix F ■ (, depend only on and hence are linear in 



spec 



(r.c) = {o,...,o,±M/ (fc) |ci; d (k) } 



(8.69) 



where V(k) are some positive constants and d(k) are the corresponding multiplicities. Here 
the index (k) labels all non-zero eigenvalues. It is put in round brackets to avoid any 
confusion with the boundary coordinate index j. Then the strong ellipticity condition 
fl3.ll ) implies 

< v (k) < 1. (8.70) 

Let us compute the fibre trace of the heat-kernel diagonal. Taking the trace of ( |6.51|) 
we get 



tT V U (t\x,x) = (4vrt)- m / 2 ic + 



c\ exp 



r{xY 



J 



I r ' 



where 



Co 



try/, d=tr v (/-2n), 



(8.71) 
(8.72) 



J(z) = tl V ${z) 



Remember that the contour C lies in the upper half plane: it comes from — oo + is, 
encircles the point u = in the clockwise direction and goes to +oo + is. 

We want to compute the asymptotics of the parametrix as r — > while t is being fixed. 
This corresponds to the limit z — ► + . The first two terms are well defined. The real 
problem is the asymptotics of the function J as z — > + . The integral over uj is calculated 
by using the formula 



duf(uj) = -27T2Res w=il/(fc) | ? | f(u) + / dwf(uj). (8.74) 

— oo 

The integrals over ( can be reduced to Gaussian integrals by lifting the denominator in 
the exponent or by using spherical coordinates. 

Note that in the integral over u all poles of the integrand lie on the imaginary axis. 
Now, if the strong ellipticity condition ( |S.70D is satisfied, i.e. all v^) < 1, then they do not 
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reach the point This is very important. This simple fact leads then to convergence of 
the integral over ( and regularity of the limit z — » + . Therefore, the heat-kernel diagonal 
is integrable near the boundary, leading to the asymptotics obtained in the previous 
Section. 

Let us instead suppose that the strong ellipticity condition is violated in that there is 
an eigenvalue v = 1 with a multiplicity d. As is shown in Ref. 1, in this case the same 
procedure leads to a singularity of the function J as z — > + , i.e. 

(m \ \ 
7 J ^' (8 ' 75) 

and hence to a singularity of the parametrix near the boundary when r — > 0, t being fixed. 
Moreover, this singularity is not integrable, which means that the L 2 trace of the heat 
kernel, Tr L 2 exp(— tF), does not exist at all! This is also reflected in the fact that the 



heat-kernel coefficients Ak/2 become singular. In other words, the standard form ( 4.17|) 
of the asymptotic expansion of the heat kernel is no longer valid. 

The singularity at the point z = results exactly from the pole at to = In the 
strongly elliptic case all poles lie on the positive imaginary line with Imu < i\(\, so that 
there is a finite gap between the pole located at the point with the largest value of the 
imaginary part and the point 

In the heat-kernel diagonal there are three types of terms now. The first class of 
terms do not vanish exponentially when r — ► oo. Those are the interior terms. They give 
the familiar interior contribution when integrated over a compact manifold. The second 
class of terms are those which are exponentially small when r — > oo and when t —>■ + . 
These are the boundary terms. When integrated over the manifold they produce the 
boundary terms in the standard heat-kernel asymptotics. In fact, these terms behave, as 
t — > + , as distributions near the boundary, so that they give well defined non- vanishing 
contributions (in form of integrals over the boundary) when integrated with a function. 
In the non-elliptic case we have however obtained also a third term. This term has an 
unusual non-integrable singularity at the boundary as r — > (on fixing t) 

ti v U {t\x,x) r ~° (4vr)- m/2 2d(m - l)r(m/2)i (8.76) 

Such a singularity is non-standard in that: i) it does not depend on t and ii) it is not 
integrable over r near the boundary, as r — * 0. This is a direct consequence of the violation 
of strong ellipticity. 

One can ask: what if the strong ellipticity condition ( |3.70| ) is violated "strongly", i.e. 
there are some eigenvalues that are larger than one, v > 1? Well, then it is not difficult 
to see that the integrals (|8.73|) defining the function J diverge for any z. Thus, in this 



case the parametrix itself, not only its functional trace, does not exist at all! 



I. G. Avramidi and G. Esposito: Heat kernel asymptotics 



16 



9 Ellipticity in Linearized Gauge Theories 

In this section we are going to show how the Gilkey-Smith boundary-value problem can 
be formulated in a general gauge theory, following Ref. 1. A linearized gauge theory is 
defined by two vector bundles, V and G, such that dimV > dimG. V is the bundle 
of gauge fields (p G C°°(V, M), and G (usually a group) is the bundle of parameters 
of gauge transformations e G C°°(G', M). Both bundles V and G are equipped with 
some positive-definite metrics here denoted by E and 7, respectively, that are Hermitian: 
E' = E,y = 7, and with the corresponding natural L 2 scalar products (, )y and (, )g- 

The gauge transformations are described by a first-order differential operator R : 
C°°(G, M) — > C°°(V, M). We restrict to the most important case when the second-order 
operator L : C°°(G, M) -> C°°(G, M) defined by L = RR, where R = 7 -1 i2 t £', is 
a Laplace-type operator with a non-degenerate leading symbol <7£,(L;£) = |£| 2 /g- This 
means that ranker^/?) = dimG. 

The dynamics of gauge fields (p G C 00 (l / , M) at the linearized level is described by a 
gauge-invariant and formally self-adjoint second-order differential operator A : C°°(V, M) - 
C°°(V, M). It is gauge-invariant in the sense that its leading symbol is degenerate and 
satisfies the identities 

a L (A)a L (R) = a L (R)a L (A) = 0. (9.77) 

We also assume that Kercr^A) = {a L (R)e \ e G G}, and hence rankcr^(A) = dimV — 
dimG. 

Instead of the operator A we introduce another formally self-adjoint second-order 
operator F : C°°(V, M) -> C°°(V : M) by 

F = A + RR. (9.78) 

Here we again restrict ourselves to the most important case when F is also a Laplace-type 
operator with a non-degenerate leading symbol o~l(F] £) = |£| 2 iy. Both restrictions made 
so far are satisfied in many interesting examples, like Yang-Mills and Einstein theories 
(for more details, see Ref. 1) 

In quantum field theory one is interested in the one-loop effective action which is 
expressed in terms of the functional determinants of the operators F and L by 

= -log Det F - log Det L. (9.79) 

On manifolds with boundary one has to impose some boundary conditions to make 
these operators self-adjoint and elliptic. In gauge theories one tries to choose the boundary 
conditions in a gauge-invariant way. Interestingly, this requirement fixes completely the 
form of the boundary operators associated to the operators F and L, respectively. 

Let us define restrictions of the leading symbols of the operators R and A to the 
boundary, i.e. 

U = a ^ A '^ B M^ v^vl&N)^, ^a L (R;0^- (9.80) 



dM' 



dM 
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Since L = RR and F = A + RR are Laplace-type operators, it follows that uu = Iq 
and II = Iy — vv. Therefore, II is a self-adjoint projector orthogonal to v, II = II, 

Tiu = m = o. 

The requirement of gauge invariance of the boundary conditions determines in an 
almost unique way that the boundary conditions for the operator L should be of Dirichlet 

type, 

e =0, (9.81) 

dM ' v ' 

and the boundary conditions for the operator F should read 

IL<p =0, Rip = 0. (9.82) 



dM 



dM 



Since the operator R in the boundary conditions (|9.82| ) is a first-order operator, the set 
of boundary conditions (|9.82| ) is equivalent to the general Gilkey-Smith scheme formulated 
in Sec. 2. Separating the normal derivative in the operator R and denoting by W$ the 
restriction of the vector bundle V to the boundary, we find exactly the Gilkey-Smith 
boundary conditions (|2.3|) with the boundary operator B of the form (|2.5| ) involving a 
first-order operator A : C 0O (Wq, dM) — > C°°(W , dM), the matrices T- 7 being of the form 

p = -vv^?v. (9.83) 

These matrices are anti-self-adjoint, P = —T\ and orthogonal to the projector II, i.e. 
IIP = PIT = 0. 

The condition of strong ellipticity then means that the matrix \(\I—iT-( = \(\I+ii>i>ij,i> 
should be positive-definite. The sufficient condition ( |3.11|) of ellipticity now reads 

Id 2 / + (r • C) 2 = ICI 2 n + (/ - n)[|c| 2 / - - n) > o. (9.84) 

Since for non-vanishing ( the part proportional to IT is positive-definite, the condition of 
strong ellipticity takes the form 

(/-n)[|c| 2 /-/i/i](/-n)>o. (9.85) 

Thus, the following theorem is found to hold: |IJ 

Theorem 2 The boundary-value problem (F, B) with the boundary operator B deter- 
mined by the boundary conditions ( p.SQ ) is gauge-invariant provided that the boundary 
operator associated to the operator L takes the Dirichlet form. Moreover, it is strongly 
elliptic with respect to the cone C — R + if and only if the matrix + ivvp,v] is positive- 
definite. A sufficient condition for that reads 

(/-n)[|c| 2 /-/^](/-n)>o. (9.86) 



We study some explicit examples of gauge theories, including Yang-Mills model and 
Einstein quantum gravity, in Ref. 1 and in another contribution to this volume. 1 13 
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